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ABSTRACT 

We consider the problem of determining the structure of the dark halo of nearby dwarf 
spheroidal galaxies (dSphs) from the spherical Jeans equations. Whether the dark halos are 
cusped or cored at the centre is an important strategic problem in modern astronomy. The ob- 
servational data comprise the line-of-sight velocity dispersion of a luminous tracer population. 
We show that when such data are analysed to find the dark matter density with the spherical 
Poisson and Jeans equations, then the generic solution is a dark halo density that is cusped like 
an isothermal (po °^ '"^)- Although milder cusps (like the Navarro-Frenk- White po °^ ') 
and even cores are possible, they are not generic. Such solutions exist only if the anisotropy 
parameter /? and the logarithmic slope of the stellar density jc satisfy the constraint jc - 2/3 
at the centre or if the radial velocity dispersion falls to zero at the centre. This surprisingly 
strong statement is really a consequence of the assumption of spherical symmetry, and the 
consequent coordinate singularity at the origin. So, for example, a dSph with an exponential 
light profile can exist in Navarro-Frenk- White halo and have a flat velocity dispersion, but 
anisotropy in general drives the dark halo solution to an isothermal cusp. The identified cusp 
or core is therefore a consequence of the assumptions (particularly of spherical symmetry and 
isotropy), and not the data. 
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. 1 INTRODUCTION 

, Hot stellar systems are held up by the stellar velocity dispersion 
' and have little or no rotation. In fact, many such stellar systems - 
giant elliptical galaxies and dwarf spheroidals - have a velocity dis- 
persion profile that is constant to a good approximation. The case 
of the dwarf sp heroidals (dSphs) has rece ived particular attention 
in recent years. iKlevna et alTfeOOll l2002h showed that the veloc- 
. ity dispersion profile of Draco is fl attish through o u t the bulk of 
the galaxy, although later work by 'Wilkinson et alj j2004) found 
evide nce for kinematicall y cold populations in the very outermost 
parts. IWalker et al.U2007h presented stellar velocity dispersion pro- 
files for seven Milky Way dSphs and found almost all to be con- 
stant to a good approximation right in to the very centre , althou gh 
the profile of Sextans seems to dip somewhat. Koch et al J ( 120073 a) 
studied Leo I and Leo II, and also found essentially flattish profiles. 

The hope has been that gathering line-of-sight velocities of 
bright giant stars in the Milky Way dSphs may provide evidence of 
the structure of the halo in these extremely dark matter dominated 
galaxies. A question of great interest is whether the dark halos of 
the dSphs are cusped, as predicted by numerical simulations in cold 
dark matter cosmogonies , or cored - but resu lts so far have been in- 
conclusive. For example. iKoch et si\ ( l2007bl) found that the veloc- 
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ity dispersion data on Leo II are consiste nt with halo dark matter 
densit ies that are both c ored and cusped. IWilkinson et al.l ^20061) 
and Gilmo re et alj ( l2007h presented an analysis of the velocity dis- 
persion data for six dSphs based on the Jeans equations and argued 
that cored density profiles were favored, partly because this also ex- 
plains the per sistence of kinemat ically cold substructure in the Ursa 
Minor dSph jKlevna et al.ll2003h and the main tenance of globular 
clusters in the Fornax dSph jOoerdt et alj200d) . 

Here, we are less concerned with modelling observational data 
on any given dSph than with understanding the generic qualities of 
the light and dark matter profiles. We consider the general problem 
of a tracer stellar population with a known line-of-sight velocity 
dispersion residing in a spherical dark matter halo of an unknown 
density law. Given the data. Sect. |2] considers what can be legit- 
imately inferred concerning the properties of the dark halo. Mo- 
tivated by the flatness of the observed dispersion profiles. Sect. [3] 
considers constant velocity dispersion solutions of the Jeans equa- 
tions. We show that the generic solution gives a halo density law 
with an isothermal cusp (po r^~)- Although other solutions are 
possible - in particular with cores or with the milder cusps pre- 
ferred by cosmologists (po r"') - they are not generic. Sect. |4] 
gives some examples, which show why Jeans modelling in the 
isotropic case has yielded results consistent with both cores and 
cusps. In our examples, however, any anisotropy drives the dark 
halo solution to the isothermal cusp. Finally, in Sect. \5\ we dis- 
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card the assumption of constancy of the velocity dispersion. Solely 
within the framework of the spherical symmetric Poisson and Jeans 
equations, we show that solutions of these equations almost always 
possess isothermal cusps, unless some very special conditions are 
satisfied either by the radial velocity dispersion or by the anisotropy 
and the logarithmic gradient of the light profile. 



2 THE JEANS DEGENERACY 

The observables are the surface brightness and the line-of-sight ve- 
locity dispersion of a stellar population. Given a mass-to-light ratio 
(Y), the surface mass density of the stellar populations (Sf ) can be 
deduced from the surface brightness. If the system is spherically 
symmetric, is then related to the three-dimensional density 

associated with the luminous material pe{r) via an Abel transform. 
Here R is the projected distance, whilst r is the three-dimensional 
distance measured from the centre of the halo. The inverse tra ns- 
form provides us with the unique pe jBinnev & Tremain3ll987h ; 
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However, even if one assumes spherical symmetry, the behaviour 
of the line-of-sight velocity dispersion does not produce a unique 
solution for the radial dependence of the radial and tangential ve- 
locity dispersions. The "luminosity weighted" (assuming a constant 
T for the stellar population) line-of-sight velocity dispersion cr^^^ is 
given by the integral 
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where j8 = 1 - (crl/crj) is the anisotropy parameter for a spher- 
ical system, and <Tr and <To are the radial and (one-dimensional) 
tangential velocity dispersions. It has long been kn own (see e.g., 
iMerrifield & Kenj|l990l : iDeionghe & Merrittlll992l) that the line- 
of-sight velocity second moment is degenerate, in that there exist 
many sets of solutions - crj(r) and /3(r) - that reproduce the observ- 
ables. 

For example, suppose that we have p[(r) from equation (|T]l. 
Then, for any given behaviour of cr^(r), the anisotropy parameter 
j8(r) can be founcQ to reproduce the observables such as 
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is the luminosity-weighted mean projected radial velocity disper- 
sion. [We have not found eq. ([3} for the anisotropy parameter in the 
existing literature.] Not all solutions are necessarily physical since 



' The first part of eq. (3j is easily verified after re-arranging eq. ^2) and ap- 
plying an inverse Abel transform. The derivation of the second part requires 
switching the order of a double integral and performing explicit integrations 
on it. 



-oo < /? < 1, which has to be checked a posteriori. However, sub- 
ject to this condition, for any given arbitrary c^ir), an anisotropy 
parameter j3(r) can be found to reproduce any observable cr^^JJR). 

Once we have peir), guessed (any) crl(r), and found the 
anisotropy parameter p{r) to reproduce the observables, then from 
the spherically-symmetric steady-state Jeans equation 

d(pecrl) 



dr 
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Here, pt is the total density such that p, = tp( + po where p( and 
Pd are the stellar and dark matter densities, respectively. With the 
self-consistency assumption, if there is no dark matter (po = 0), 
and the mass-to-light ratio T is additionally constant, then the 
choice of al{r) is uniquely determined by the coupled Poisson 
and Jeans e quations (s ee e .g., iBinnev & Mamonlll982l : lTonrvll983l; 
IDeionghe & Merritt 1993). 

However, if T is varying or p^ir) + 0, there is no unique 
choice of cr\(r). That is, any a\{r) is allowed subject only to the 
constraints that -oo < p(r) < 1 and pt(f) > 0. Consequently, with- 
out further observational constraints or simplifying assumptions, 
determining a model to reproduce the observed Sf(/?) and a\^JJi) 
is completely indeterminate in the spherical case. 



3 JEANS MODELLING WITH CONSTANT VELOCITY 
DISPERSIONS 

3.1 Isotropy 

The simplest assumption to make is that of isotropy {fi = 0). Then 
crl(r) is recovered uniquely from o-^^^(R) by an inverse Abel Trans- 
form. The Jeans equation reduces basically to hydrostatic equilib- 
rium with the "pressure" being equal to P = pecr^, that is. 



V(pfO--) = -piVi/i 



(5) 



where cr is the one-dimensional velocity dispersion of the tracer 
population and ifi is the gravitational potential. 

If we further assume that o-^^^(R) is a constant cr^, as suggested 
by most of the observations, then the unique solution is cr^(r) = 
(Tq. Under this assumption, the central properties of the halo are 
severely restricted, as we now show. Since V(pfO"-) = cr^Vpc, the 
Jeans equation indicates that Vpf and Vi^ are (anti-)parallel every- 
where. This further implies that the surfaces of constant pf and ij/ 
coincide, and thus p( can be considered as a function of i/^. Conse- 
quently, Vpf = {dp(/difi)Vi// and equation ^ reduces to 

By solving this differential equation, we find that 



(6) 



Pe = Po exp 
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where iJ/q = cTq Inpo is an integration constant. Combined with the 
Poisson equation under the assumption that the potential is gener- 
ated by the dark matter halo of a density po (i.e., pn, » pf and so 
Pt ~ Pd), we obtain 
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This is an interesting equation, as the dark matter density, which we 
wish to know, depends on the Laplacian of the luminosity density. 
The combined integro-differential equations l[T]l and ([S} now relate 
the dark matter density to the observables, ctq and Sf (-R). 
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The implications of equations ([7} and ([8} on the behaviour of 
tracer populations in dark matter halos are of considerable interest, 
even before we apply to any specific example. First, equation ([7J in- 
dicates that the potential is finite for any finite luminosity density. 
Consequently, we find that any cored luminosity density implies 
that the central potential well cannot be infinitely deep and so the 
dark matter halo density profile is also cored, or diverges strictly 
slower than the singular isothermal sphere (r^^) if it is cusped. 
Similar argument also leads us to the conclusion that any cusped 
luminosity density would only be supported by a cusped halo di- 
verging at least as fast as a singular isothermal sphere. In fact, this 
latter conclusion can be sharpened with further analysis. Under the 
assumption of spherical symmetry, equation ^ may be written to 
be 

dye AnG ^ 

7c + r-^ = —rPor (9) 
dr 

where yc = -(dlnpf/dlnr) is the logarithmic slope of the density 
of the tracers. By taking the limit to the centre (r — » 0), we find that 

limPD'-' = /^7f.o, (10) 

where yf o is the limiting value of yr towards r — > 0, that is, the log- 
arithmic cusp slope of the luminous tracer density. This indicates 
that if the luminosity density is cusped with yf o > 0, the halo den- 
sity must be cusped as po ~ ''"^ like a singular isothermal sphere. 
On the other hand, any cored p[ indicates that y^-.o = limr_,o(pD''^) = 
and thus the halo density may not diverge as fast as or faster than 
the cusp of the singular isothermal sphere. Finally, pf with a central 
hole implies th at yt^p < and is th erefore unphysical as it will lead 
to Pd < (c.f. lAn & Evanil2006h 

3.2 Anisotropy 

In reality, the velocity dispersion tensor of a "coUisionless" stellar 
system is not necessarily isotropic. However, provided the gravi- 
tational potential is still spherically symmetric, the Jeans equation 
reduces to 

Id/ ,x dl/r 

where / = exp J(2y6/r)dr is the integrating factor (e.g., / 
is constant). 

Inspired by the discussion in the preceding section, we now 
consider the case of constant radial velocity dispersions with an 
arbitrary functional form of fi. Note that it is possible, from equa- 
tion (O, to find y8(r) that is consistent with the observed Ef(i?) and 
(T^g^iR) once we ascribe a particular behaviour to cr-^(r). Hence, we 
can in principle find such a model that produces a constant line- 
of-sight velocity dispersion (or any other desired form), which is 
indicated by the observations using equation (O. 

If (T^ = CTp is a constant, equation dl It becomes 

dif/ T d 
d7 = "'^"d'r 

and consequently we find that 
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Here, the last equation can be re-cast similarly to equation (O so 
that 
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This is basically the same as equation l[9j, except yc is replaced by 
7f - 2/3. We note that the result does not require the assumption 
that yS is a constant. The implication of this result is quite similar 
to that of equation At the limit towards the centre (r — > 0), 
we find that yi-g > 2/3o implies that pu ~ r^^ whilst we infer 
that limr_,o(pD''^) = if yc,o = 2/3o. Here, is the limiting valu e 
of anisotropy parameter at the centre. As per I An & Evans! ( l2006h . 
7f.o < 2y8o is unphysical (although pf o does not self-consistently 
generate iff, the potential well depth is finite at the centre so that 
their result holds) since it indicates negative halo density. 

In summary, therefore, given a tracer population with a con- 
stant radial velocity dispersion, the generic solution of the Jeans 
equation for the dark matter is cusped like a singular isothermal 
sphere (po oc r^-}. Milder cusps (like pn r^' ) and cores are 
possible, but they are not generic. Such solutions only exist if the 
anisotropy parameter /3 and the logarithmic slope of the stellar den- 
sity ye satisfy y^.o = 2y3o. 



4 EXAMPLES 

4.1 A Plummer Light Profile 

Plummer's law is commonly used to model the light of dSphs (e.g., 
lLakelll990l : IWilkinson et"al]|2002h . Assuming a constant mass-to- 
light ratio Y, then the surface density is 



So 



(15) 



{\+R^lrl)^ 

Here, is the radius of the cylinder that encloses half the light, 
whilst the total luminosity is L = Ttr^So/Y. It is straightforward to 
establish via equation l[T) that the stellar density is 

320 



1 



(16) 



(17) 



4ro (I + r2/r2)5/2 
Now, using equation ([8j, the dark matter density must be 

47TGr2 (I +r2/r2)2' 

which is a cored isothermal sphere (see lEvanslll993l) . However, as 
the model is isotropic (j8 = 0) and the stellar density is cored at the 
centre (y^ o = 0), this dark halo solution corresponds to the special 
case yf o = 2^60. 

Now suppose the assumption of isotropy is dropped. Using 
equation we find that the dark halo density acquires an addi- 
tional term that behaves near the centre as 



(18) 



2nGr^ ' 

from which we deduce that ySo ^ so that the model is tangentially 
anisotropic as r — > 0. Notice that this has changed the behaviour of 
the density at the origin - the halo law has now an isothermal cusp. 
This is in accord with the general result that provided yt^ > Pq, the 
cusp is isothermal. 

4.2 Exponential Light Profiles 

Another set of profiles often used to model dSph light distributions 
is based on the exponential l aw ( ISersic|[T968l : iFaber & Linlll983l : 
llrwin & Hatzidimitriou|[l99^ . Suppose that the three-dimensional 
density law pe{r) is exponential 
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M'-)=Po exp|---j . (19) 

The corresponding surface density is cored as 

EKW) = Eo (20) 

where Eq = 2rdPo and Ky{x) is the modified Bessel function of the 
second kind and of order v. We refer to this as the Bessel profile. 
The half-light radius is ~ l.Qllri, whilst the total luminosity is 
L = 4T:rjl,o/T. The dark matter density inferred from equation {8} 
is then 



Pd 



(21) 



which is the r"' cusp beloved of cosmologists (e.g. jNavarro et al] 
Il996h . In fact, a three-dimensional density law pi = poe leads 
to an infinite dark matter cusp of form pn oc r"'-""'. 

Suppose instead the surface brightness profile is modelled 
with an exponential law 



I.((R) = So exp - — 



(22) 



The luminosity density is then (see e.g., lKenj|l992h 

which is logarithmically divergent at the centre. The dark matter 
density is now 



(23) 
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K„{r/R^) and n 
1 



: 0, 1 . It is singular as r — » 0, 
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(25) 



which exhibits strictly slower divergence than a singular isother- 
mal sphere. In other words, once the luminosity density has been 
assumed to be of exponential form (either in projection or in three- 
dimensions), then the inferred dark matter density is cusped, but 
the cusp is always weaker than the isothermal cusp po ~ r^^. 

It is easy to see that equation (|9} still applies since — > 
as r — » 0. Again, the isotropic models are somewhat unusual - 
the introduction of anisotropy drives the dark halo solution towards 
an isothermal cusp. The fact that the terms in the density and the 
anisotropy decouple in equation l ll3t means that the assumption of 
any central anisotropy gives the same additional contribution to the 
halo density l |18l l as in the previous example. 



5 THE GENERAL CASE 

Let us now gain insight into the general case by discarding the 
assumption that the radial velocity dispersion is constant. We de- 
rive the extension of our result to the general spherically symmetric 
case. By rewriting the spherical steady-state Jeans equation l|4j, we 
obtain (under the assumption that p, = po) 

f 9 ,/ diner;'! 

47tG J po(f)r^df = rcrjiy(-2/i- 



and 



AnGp^r- = cr^ 1 + 
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Now note that if the dark matter density has an isothermal or 
steeper cusp, then the left-hand side of equation ill) tends to a non- 
zero value as r — > 0. However, if the dark matter density is cored 
or diverges more slowly than r^^, then the left-hand side vanishes 
as r — » 0. Then, for the right-hand side also to vanish as r — > 0, one 
of the following three conditions must hold at the same limit 



(i) 



din (tJ 
dlnr 
dlno-? 



din r 
(iii) 0-2 ^ 0. 



-1 



rr - 2/3 



Case (i) implies that crj ~ r"' . Excluding the possibility that there is 
a black hole at the centre, then the central potential must be finite or 
divergent strictly slower than the logarithm (as poi^ — > 0), and the 
velocity dispersion diverging as a power law cannot be supported. 
Case (ii) implies that crj ~ r'"^^^. If yc > 2fS, then o-J — » and so 
this may be subsumed into Case (iii). If yc < 2/3, then cr^. diverges 
as a power law and so is again unphysical. This leaves only y/ = 2(3 
as an independent possibility. 

Consequently, we have established that if the dark matter den- 
sity is cored or falls off with a cusp less severe than r"^, then either 
y( = 2y3 or 0-? = at the centre. The discarding of the assumption of 
constancy of the radial velocity dispersion permits one additional 
possibility, namely that o-? falls to zero at r = 0. 

In summary, therefore, we have a surprisingly strong and gen- 
eral result. Given a tracer population in the spherical Jeans equa- 
tion, the generic solution for the dark matter is cusped like a sin- 
gular isothermal sphere fpo oc r^^). Milder cusps and cores are 
possible, but they are not generic. Such solutions exist either if the 
anisotropy parameter p and the logarithmic slope of the stellar den- 
sity yt satisfy yi^ - 2/3o or if the central radial velocity dispersion 
cr^^Q vanishes. 

Such a strong result may seem astonishing. It is worth remark- 
ing that this result is essentially due to the spherical symmetry as- 
sumption, and in particular, the coordinate singularity at the origin. 
A real dSph is not likely to be exactly spherically symmetric at the 
centre. So, it may be possible for a dSph to have a finite velocity 
dispersion towards the centre by mild symmetry-breaking. Outside 
of the very central region, a dSph can often be well-approximated 
by the idealized spherically symmetric fiction. Our theorem ther- 
fore really lays bare the dangers of inferring the central density 
profile from the almost universally made assumption of spherical 
symmetry. 



6 SUMMARY 

We have studied the problem of deducing the dark halo density 
from the surface brightness and velocity dispersion profiles of a 
tracer population. This has immediate application to the dwarf 
spheroidals satellite galaxies (dSphs) of the Milky Way. 

If the stellar population generates the gravity field, then the 
Jeans and Poisson equations give a unique solution for the density, 
the mass-to-light ratio and the anisotropy of the spherical model, 
consistent with the observed surface brightness and the line-of- 
sight velocity dispersion (Binney & Mamon 1982; Tonrv 1983). Of 
course, this does not apply to the case of dSphs, in which the den- 
sity of the stellar population is dominated by the dark halo density. 
Now, the problem suffers from the well-known mass-anisotropy de- 
generacy. 

The line-of-sight velocity dispersion profiles of the Milky Way 
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dSphs appear to be usually flat (see e. g., lKlevnaetal.ll200ll , l2002l : 
iKoch et al.ll2007al . IWalker et id]|2007h . which suggests the simple 
assumption that the velocity dispersion tensor is isotropic and has 
a constant value. Then, any inference as to the central behaviour of 
the dark matter potential is controlled by the assumption as to the 
light profile of the tracer population. If the light profile is cored, 
then a dark matter halo density that is itself cored is deduced from 
Jeans modelling. If the light profile is cusped like an exponential 
law (or its variants), then a dark halo density that has a milder 
cusp than isothermal (such as the Navarro-Frenk- White cusp of 
Pd cx r"') is deduce d. This provides the explanation as to why pre- 
vious investigators ( Koch et al ]|2007bl J Wilkinso n et al. 2006) have 
concluded that the data are consistent with both cusps and cores. 

However, velocity anisotropy in the stellar population has a 
dramatic effect on the dark halo density recovered from Jeans mod- 
elling. If a tracer population has a constant radial velocity disper- 
sion, then the generic solution for the dark halo is always cusped 
like a singular isothermal sphere (pn oc r"-). Milder dark mat- 
ter cusps (like pu oc r"') and cores are possible, but they are not 
generic. They can occur only when the condition yc = 2/3 is ful- 
filled at the centre, where /? is the anisotropy parameter and jf is the 
logarithmic slope of the stellar density. Note that many of the com- 
monly used dSph models (such as Plummer or exponential profiles 
with isotropic velocities) correspond to the special case yf = 2/3 
and so any conclusions inferred as to the dark halo law may not be 
beyond reproach. 

Finally, even if the assumption as to the constancy of the ra- 
dial velocity dispersion is discarded, then almost the same theorem 
holds true. If y( = 2p at the centre or if crj falls to zero at the centre, 
then dark matter cores and milder cusps than isothermal are pos- 
sible. The generic solution, however, remains the isothermal dark 
matter cusp, at least within the framework of the spherical symmet- 
ric Jeans and Poisson equations. 

Here, our examples and analysis have shown how Jeans solu- 
tions may be telling modellers more about their assumptions rather 
than the theoretical implications of the data. Of course, there is in 
principle more information in the discrete velocities than in the ve- 
locity dispersion profile and the Jeans equations. The best response 
to the degeneracy of the problem is to seek further observational 
constraints (perhaps h igher moments from the line profile, see e.g., 
iMerrifield & Kenllll99Q) or additional insights on the behaviour of 
the anisotropy (maybe from simulations or from a detailed analysis 
on the underlying physics). 
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